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Abstract 

oo ■ 

The laser-ion acceleration in the ultra-short and ultra-intense laser-matter interactions attracts more and 
more interest nowadays. Since electrons gain relativistic energy from laser pulse in a period of several 
femtoseconds and driven away by the ponderomotive force of laser pulse, a huge charge-separation field 

S: 

pulse is generated. In general cases, the ion acceleration is determined by this charge-separation field. A 

c3 , 

O ,| novel general time-dependent solution for laser-plasma isothermal expansions into a vacuum with different 

types of the scale length of the density gradient which correspond to different charge separation forms is 
obtained. The previous solutions are some special cases of our general solution. A series of new solutions 
have been proposed and may be used to predict new mechanisms of ion acceleration. However, many 
^ ' unaccounted idiographic solutions that may be used to reveal new acceleration mode of ions such as shock 

^ . wave acceleration, may be deduced from our general solutions. 
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I. INTRODUCTION 



The generation of energetic proton and acceleration mechanisms in the ultra-intense laser 
pulses interaction with thin targets attract more and more interest nowadays Ilil2il3i|4j,l5|]. Their 
progress can provide fundamental theory for inertial confined fusion (ICF) and promote the real- 
ization of it effectively. The ultra-short and energetic ion beam allows for an increase of energy 
resolution in the Time Of Flight experiments, the investigation of the dynamics of nuclear pro- 
cesses with high temporal resolution and the study of spallation-related physicsflg]. 

When a relativistic laser pulse interacts with a plasma, the laser-produced fast electrons with a 
unique temperature, k B T e , determined by the laser ponderomotive potential are instantly created 
and then driven away. However, the ions are still resting due to the large mass and then a high 
charge separation field generates. Furthermore, the plasma is assumed isothermal since the con- 
tinuous energy supply of the laser pulse in the pulse duration. No matter proton shock acceleration 
(PSA)|Q] in laser-plasma interactions or target normal sheath acceleration (TNSA)[4, 5, 8, ^] and 
so on, the ions are accelerated by high charge- separation field. The key point of ion acceleration 
is the hot-electron density distribution which decides the spatial and temporal distribution of the 
charge- separation field. 

In this paper, a general solution for plasma isothermal expansions into a vacuum is proposed 
with the assumption: the ion density distribution can be represented a function with separable 
variables in the transformation system. With the solution, the separate charge distribution, electric 
field, electron velocity, ion velocity, and fronts of ions and electrons are all predicted. For different 
scaling length of the density, the solution corresponds to different expansion mode of plasmas. In 



nn 



some special cases, the solutions have been achieved by previous pursuers (5|, IE 1011 . A series of 
new special solutions have been described and the corresponding acceleration modes have been 
discussed in detail. It is pointed out the shock wave forms for some types of the plasma density 
gradient and large scale length. 



H. TIME-DEPENDENT ION ACCELERATION DUE TO STRONG CHARGE SEPARATION 

For convenience, the physical parameters: the time, t, the ion position, /, the ion velocity, v, 
the electron field, E, the electric potential, </>, the plasma density, n, and the light speed, c, are 
normalized as follows: i = u) p ,ot, I = I /Aw, u = v/c s , E = E/E , (p = e(p/k B T e , h = n/n e0 , c = c/c s , 
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FIG. 1: (Color online) The ion velocity, and electron velocity, n e VS the self-similar variable, for 
L = -f}\ = -1 and N\ = 1 in two cases: u$ - c 2 - 0.02 ; uq - c 2 - -0.02. 

where n represents n t (or n e ) which is the ion (or electron) density, n e0 is the reference hot-electron 
density, c s = ^Zk B TJmj is the ion acoustic speed, co piQ = ■^Zn e oe 2 /m i eo is the initial ion plasma 
frequency, Ado = c s <x> pi o, c is the light speed and E = k B T e /eA D o- Here e is the elemental charge. 

The reference frame used here is r = t, g = x/t. With the transformation, the equations of con- 
tinuity and motion are obtained easily in the new coordinate system. The ion density is assumed to 
satisfy: n ; (r, , £) = Ni(t)N 2 (0- Here, N\ = 1 corresponds the self-similar ion density and solution. 
If N 2 = exp(-g/f3\ - 1), where fii is a constant, the self-similar solution is for a neutral-plasma 
isothermal expansion into a vacuum given by Huang et al. |5|] for the impurity ions. /3i = 1 
corresponds the classic self-similar solution given by Mora in [|80. If N 2 = n (^/^o)~ 2/a , where 
a e (0, 2), the solution is for a half-self-similar non-neutral plasma isothermal expansion into 
a vacuum proposed by Huang et al. Jlp| . For dNi/dr ± 0, the analytic solution has not been 
reported. That is what given out by us next. 

Combining the continuity and motion equation of ions gives the general solution of the ion 
velocity and potential in the ion region: 



Ui = Of 2 (£) - tfi(T)0fi(£), 

(p = -a£/2 + /[<*» + (Si + S 2 )ai]d^, 



(1) 



where Si = t/L t , S 2 = r 2 /L Ti2 , L T = [dln/Vi/dr] 1 is the time scale length of the ion density, 
Ni(t), L Ti2 = [d^nNJdT 2 ]- 1 , a x = F f F~ l dg and a 2 = F f F~ 1 ^' '/Ldg ', a 3 = 5 Y ai + £ - a 2 , 
F = exp"i^ /i and L = [d\nN 2 /dg]~ l is the scale length of the time-independent ion density, 
N 2 (ij). And the electric field is E = -(Si + S^ajr - o|/Lt + Sia^/r. 



Combining Eq. (OQ) and Poisson's equation, the electron density satisfies n e = rii - 8n, where 6n 
is decided by: 

2of(l+ L'/2) (2 + 35 1 )a 3 -(5i+5 2 )a 1 -((52-5 2 )L 
* = L? • (2) 

where L = dL/d£. From Eq. ©, lim,-^ £n — » and it means that the plasma tends to neutral as 

r — » +00, whatever the initial state is. 

We will confirm the ion front and electron front with physical discussions and on the bases of 

Poisson's equation and the continuity of the potential and electric field the next. The first case is 

the electron front is beyond the ion front. Therefore, beyond the ion front, the ion density is zero 

and it can be assumed that the expression of electron density is a smooth expansion with respect 

to self-similar variable, With the expression of electron density, solving the continuity equation 

of electrons, the electron density is u e = U{ + 6u, where 8u satisfies: 

-6na 3 - (da 4 /dT)/T + a A /r 2 
ou = (3) 

n e 

where a 4 = (61 + 5 2 )a\ + a 2 /L - 6\a 3 and then Tda 4 /r = (61 - 6 2 - S3 + S\ + 6 l 6 2 )ai + 2a 3 (£ - 
a 2 - 8 2 ai)/L + (S 2 - 5i)a 3 , and 6 3 = T 2 d 3 NJdr 3 . In fact: 6n = -(da 4 /d£)/T 2 . With Eq. ©, 
lim^+oo u e -> U{. 

Therefore, with the electron density and Poisson's equation, the electric field beyond the ion 
front satisfies: 

E(t,$ = -tN 1 I N 2 df-^, (4) 

where represents the value of £ at the ion front. In the ion region, the electric field is E(t, = 
-a 4 (^)/r. FromEq. ©, satisfies: 

Ni N 2 d£ = f-, (5) 



where % e j stands for the position of the electron front there the electron density is zero. 
With the expression of the ion front, maximum ion velocity is given by: 

Ui, m = a 2 {^ u \ Tacc ) - 5iari(£,/| Tacc ), (6) 

where r acc is the acceleration time, which is about 1-2 times of the laser pulse duration for the 
ion acceleration in the laser- solid interactions. 

If the electron front is before the ion front, the ion velocity is larger than that of electrons. 
However, in reality, this situation can not happen. Therefore, we ignore the solutions in this case. 



Before we discuss the different cases of the scale length, L, a very special solution that does 
not rely on the choose of L is given. If N\ oc 1 /r, the special solution is: 

m = n e = Ui = u e = ft (p = (f> , E = 0, (7) 

This solution corresponds a neutral plasma with no charge separation and a constant velocity 
expands into a vacuum. 

In two special cases, the special solutions are familiar to us. 

Case one: the scale length of ion density, L, is 0-degree polynomial in ft L = -fix. If L = -fix, 
F = exp^' +1 , N 2 = F 1 and u = u{£ = ft) at ft = -fii- Then a v = -0^1-aF), a 2 = £+fix +c 2 F, 
a 3 = +(5i(l -CiF)] -c 2 F and u t = + c 2 F + J S 1 (5 1 (1 -ciF), where Ci and c 2 are integral 
constants. 

For c\ = 1 and c 2 = 0, uq = 0. If 5i < exp~^ //3l_1 , ions are accelerated. Oppositely, if 8\ > 



exp ^ l , ions are decelerated. Therefore, as Huang et al. (4j, LLJJ] pointed out <5i « 1 > exp 



-f/A-l 



for ^ > ft = -j8i , this solution is not suitable to describe the ion acceleration in the ultra-intense 
laser-foil interactions. 

For c\ = c 2 = 0, a special solution: a\ = -fix, a 2 = ^ + fix, a 3 = -fii[l + 6x], u t = ^ + fix + fixSx 
and 5n = 0. Since = -fi\(\ + 5 2 + 26 1), assuming ft,/ = +oo, from Eq. ©, the ion front is 

tij=PM - - )-fii. (8) 

1 + 26 1 + 6 2 



The main part of Eq. © is the similar as that given by Huang and co-workers lllin using physical 
discussion. This solution is a special time-dependent solution for neutral-plasma isothermal ex- 
pansion into a vacuum given by Huang et al. Ill ill . In special case: iVi = 1, with Eq. ©, the ion 
front is governed by: ft/ = fix[\n(r 2 ) - 1]. Considering the initial conditions, the results given here 
are the same as that given by Huang et al.[5] and Mora Q8D (where fix = 1). Huang et al. and Mora 
obtained the results through the physical discussion about the Debye length of electrons instead 



of analytic deductions. However, the same analytic method has been used by Huang et a/ J10h to 
deduce a special result of the following solution in case two. 

For Nx = 1, Uj = ij + fix + c 2 F. With Eq. © and Eq. ©, 6n = -2c 2 F(fi x + c 2 F)/r 2 and 
6u = [(j3x+c 2 F) 2 (2c 2 F + \/fix)]/[T 2 F- l +2c 2 F(fix+c 2 F)]. For« = c 2 = (ft = -fix), 6n = and 
the solution corresponds the general self-similar solution for neutral-plasma isothermal expansions 

u 

into a vacuum pointed out by Huang et al. [5] . However, if w = c 2 > (ft = -fix ), the acceleration 
described by the solution , which is not for a neutral-plasma expansion, is more efficient than the 



classic solution. If uq = c 2 < -j8i (£o = -fii), the ion will be accelerated in the opposite direction. 
If -fii < uq = c 2 < (^o = -fii), the ion will be accelerated first and decelerated then. The density 
distributions for c 2 > and c 2 < have been shown by Fig. [TJ 

Fig. Q] shows: u = c 2 = -0.02 and the electron density is smaller than the ion density for £ < 3 
and the ion acceleration is less efficient than that in the neutral-plasma case; u = c 2 = 0.02 and 
the electron density is larger than the ion density for all £ and the ion acceleration is more efficient 
than that in the neutral-plasma case, since the hot-electron number is so much more than that of 
ions. In these two cases, the density difference becomes small with the increase of r and tends to 
zero as r — » +oo. The electron front is determined by: u e (j; e j) « c in the relativistic laser intensity. 
Then, with g e j, the ion front can be obtained by Eq. ©. 

For c 2 = 0, it is equivalent to: g e j+/3i +/S i exp(^ e j//3 i + I) It 2 ] « c. With this, the electron front 
satisfies: {j e j < c -/3i(l + exp(l)/r 2 ). However, in previous works given by Mora QSD and Huang 
and co-workers [5], the electron front is taken as the positive infinity. Then, with and Eq. ©, 
the ion front can be obtained by: = /3iln[r 2 /(l + jS'j] - fii, where X = fii/(c - g e j - fix). 

Case two: L is an one-degree polynomial in L = -atj/2 and a e (0, 2). In this case, 
F = |^ | 2/ff = ? la , a, = -P£ Q Q + ciF), a 2 = (1 + + c 2 F), a 3 /£ = + 50 - 

[6\Pc\ + (1 + /3)c 2 ]F, where /3 = a/(2 - a), C\ and c 2 are integral constants. With Eq. (QQ), 
Ui/^Q = £[1 + J3(l +6i)] + [Sifici + (1 + /3)c 2 ]F. With this, the ion acceleration is determined by the 
initial state and the density distribution (neutral or not). 

For u Q {% = &) = 0, a = c 2 = -1, and the ion velocity is + P(l + - < for 

£ > > 0, u, ■ = for £ = £o , and u, > for < £ < £o- All the ions move to the central point: 
£ = and the farther the distant, the larger the ion speed. With Eq. ©, the density difference, Sn, 
is calculated for N\ = 1. In this case, the electron front is given by 6n(£ eJ ) = f eJ 2/a . WithEq. ©, 
the ion front is obtained. Fig. |2fa) shows the ion and electron density at r = 1 and a = 0.5, 1, 1.5. 
The ion and electron front are all obtained. This situation is not an efficient acceleration mode 
since maximum ion velocity for any a e (0, 2) is very finite as shown by Fig. Ob). 

However, if c 2 = 0.02 and Ni = 1, the ion velocity is + 0(1 + + 0.02^). Fig. g] 

shows the ion acceleration is efficient in this case, especially for a = 1.5 and £o = 1- The electron 
and ion front are obtained with n e (^ e J) = and Eq. © and shown by Fig. ©(a). 

For ci = c 2 = 0, ut = £ + #(1 + tf^and = -(1 + /?)/^ 2 (l + 8 i ) 2 /2 - (6 2 - 5 2 )/^ 2 /2. 
5n = [fi8 2 + /3 2 6 2 + 0(1 + P){28\ + l)]/r 2 This solution is also time-dependent and can be used to 
describe the energetic ion acceleration with an enhanced electron tail. For iVi = 1, the solution has 
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FIG. 2: (Color online) (a) The electron and ion density VS the self-similar variable, for L = -(a/2)g 
in three cases: a = 0.5, 1, 1.5. The electron fronts are all before the ion fronts, (b) The ion velocity VS 
the self-similar variable, for L = -(ar/2)£ and a = 0.5, 1, 1.5. In Fig. |2ta) and (b), uq(^ = £o) - 0, 
ci=c 2 = = = I- 
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FIG. 3: (Color online) (a) The ion and electron front VS the expanding time, t; (b) the ion velocity VS 
the self-similar variable, for L = -fig/2 and N\ - 1. Here, C2 - 0.02, a = 1.5 and £o - 1- The ion is 
accelerated efficiently and about 6 times acoustic velocity at r = 10. 
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been obtained by Huang et al. llOD. a 4 = -[6 2 + S { ) 2 + 1 + 26^ and £•/ = (a/lf&j, and 
the electron front is given by: 

geJ { [6 2 +/3(l+ d l ) 2 + l + 26 l ]/3 1 ' K) 

In special case: iVi = 1, the results from Eq. © are the same as that obtained in 

Two special cases: L are zero-degree and one-degree polynomials in £ have been calculated 
and some time-independent solutions are the same as the previous works. For different L, the ion 
accelerations are different. The key variable is the scale length of the time-independent electron 
density, L. In the neutral-plasma case, L = -fii (J3\ e (0, 1)), which is a zero-degree polynomial in 

In the hot-electron-tail case, L = -(a/2)g (a/2 e (0, 1)), which is a one-degree polynomial in 
Even for the same L, different integral constants induce different solutions. However, the essential 
determinant is the charge separation: 6n. Therefore, it is concluded that the charge separation of a 
plasma determines the ion acceleration. 

Here, the case: L is a quadratic polynomial in -/3 2 £ 2 (J3 2 e (0, 1]), is considered the first time. 
For a n-degree polynomial: L = can be considered in the same way. 

If L = -J3 2 e(j3 2 6 (0, 1]) and N x = 1, F = exp(-l/M + 1/Mo), *i = F~ l and u t = a 2 = 
F f F~ l d£ +u F, where u = = £ ) is a constant. With this, « 3 = £-a 2 , a 4 = -(1 -a 2 /^) 2 //3 2 , 
5n = 2(1 - « 2 /£)(l + fi 2 a 2 - a 2 l^)l0\T 2 . Since lim^ +00 a 2 —* jj + u exp(l /fotjo), some results are 
obtained: lim f ^ +00 ut — » £ + woexp(l//3 2 £o), lim^ +cx > ffs — > 0, lirn^+oo qt 4 — » and lim f _> +0o 5n — > 
0, lim^ +0o rc ; — » n e . With Eq. ©, the ion front and the electron front are all positive infinity. 
However, dui/dg < for any u and £ large enough. Therefore, the shock wave forms for £ large 
enough in this case. 

Similar with above discussion, for L = -J3 n ^ n , the physics of the expansions of plasmas can 
be obtained easily and the same results of them are: the plasma front is positive infinity and the 
boundary condition should be added exteriorly. In fact, it is easy to prove that the analytic form of 
the solutions can not be obtained and they may correspond to shock section for L = n > 2 

and £ large enough. 

It is required that dL/d£ < and L + n > 2 in order to make the limitation of the plasma 

density be finite as £ tends to positive infinity. L = -f3\ and L = -atj/2 are two simple cases. 

In order to show time-dependent solutions, it ia assumed that iVi = kt oc t, t e [0, t u ] as an 
example in two cases: L = -fa and L = -ag/2. N\ = kt, then 8\ = 1 and 5 2 = 0. For L = -f}\ 
and c\ = c 2 = 0, u t = £ + 2/3 { and w,-,/ = ySiln(/cr 3 e/3), where e = 2.718.... For L = -ag/2 and 
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Cl = c 2 = 0, Ui = (1 + 2/3)£ and % eJ = [Kn Q /(4j3 + 3)(3] a/2 T 3a/2 . Considering the dependence of the 
electron density on r described by Huang et aim [4], these solutions can be used to describe the 
influence of the hot-electron recirculation on the ion acceleration at the rear of the target heated 
by ultra-intense laser pulse. Similar with above discussions, the ion front and electron front in the 
general case can also be obtained. Here, it is not repeated again. 



HI. CONCLUSION 



In conclusion, a general time-dependent isothermal expansion for the ion acceleration due to 
charge separation is proposed on the base of the equations of continuity and motion of ions and 
Poisson's Equation. As examples, several new solutions for each types of the plasma density 
gradient have been proposed. Especially, for L = n > 2 and £ large enough, we pointed out 

that the shock wave solution exists. However, the analytic formation can not be achieved easily 
here. 
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